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Here we would like to review recent progress on the resummation on logarithms of Os in heavy quarkonium. We 
will mainly focus on the phenomenological relevance of these achievements. Determinations of the »7i,(15') mass, 
Bc{lS) hyperfine splitting, inclusive electromagnetic decays and implications for t-t production near threshold. 



1. Introduction 

Heavy quark-antiquark systems near threshold 
are characterized by the small relative velocity 
V of the heavy quarks in their center of mass 
frame. This small parameter produces a hier- 
archy of widely separated scales: m (hard), mv 
(soft), mv'^ (ultrasoft), ... The factorization be- 
tween them is efficiently achieved by using effec- 
tive field theories, where one can organize the 
calculation as various perturbative expansions on 
the ratio of the different scales effectively produc- 
ing an expansion in v. The terms in these series 
get multiplied by parametrically large logs: Inw, 
which can also be understood as the ratio of the 
different scales appearing in the physical system. 
Again, effective field theories are very efficient in 
the resummation of these large logs once a non- 
relativistic renormalization group (NRG) analysis 
of them has been performed. We will review in 
this paper recent progress on the resummation of 
the above logarithms, within the context of pN- 
RQCD [1], in the weak coupling regime {v ^ as). 
Besides the pure theoretical interest of these com- 
putations, they may also have an important phe- 
nomenological impact in several situations. Let 
us enumerate a few of them. The determination 
of the bottom and charm masses (using the exper- 
imental value of the ground state heavy quarko- 
nium masses or non-relativistis sum rules). The 
determination of the ?7f,(lS') mass, the hyperfine 
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splitting (HFS) of the ground state Be system, 
or theoretical improved determinations of the r]c. 
One can also try to obtain improved determina- 
tions for the inclusive electromagnetic decays of 
the heavy quarkonium. On the other hand the ap- 
plication of this program to t-t production near 
threshold at the Next Linear Collider is one of 
the main motivations to undergo these computa- 
tions. In this paper we review the phenomeno- 
logical analysis already available in the literature 
and outline possible future work. 

2. Hyperfine splitting 

Analytical expressions for the HFS of heavy 
quarkonium (for the equal and non-equal mass 
case) are available with leading log (LL) [2,3] and 
next-to-leading log (NLL) [4,5] accuracy. 

For the case of bottomonium, these results have 
been used in Ref. [4] to give predictions for the 
mass of the rib(lS) (using the experimentally very 
well known value of Mx(is)) with great precision 

M{r]b{lS)) = 9421 ±ll(th)t^((5a,) MeV, (1) 

where the errors due to the high-order perturba- 
tive corrections and the nonperturbative effects 
are added up in quadrature in "th" , whereas 
"^tts" stands for the uncertainty in as{Mz) — 
0.118 ± 0.003. This prediction is of great experi- 
mental interest. The discovery of the r]^ meson is 
one of the primary goals of the CLEO-c research 
program [6] and there are experimental proposals 
for its detection at Tevatron too [7]. It has also 
been argued that the HFS can be used to search 



for new physics [8]. Therefore, an accurate pre- 
diction of its mass M{rjb) is thus a big chaUenge 
and a test for the QCD theory of heavy quarko- 
nium. For instance, this prediction can be com- 
pared with those obtained either in lattice [9] , po- 
tential models (see for instance Ref. [10]) or sum 
rules [11]. It seems to be a general trend that our 
result is larger than the lattice predictions and 
smaller than most of the potential model results. 
We would also like to remark that the inclusion of 
resummation of logarithms has a sizable effect in 
the determination of the rii,{lS) mass. We illus- 
trate this fact in Fig. 1 from Ref. [4] . In this fig- 
ure, the HFS for the bottomonium ground state 
is plotted as a function of /i in the LO, NLO, 
LL, and NLL approximations. As we see, the 
LL curve shows a weaker scale dependence com- 
pared to the LO one. The scale dependence of 
the NLO and NLL expressions is further reduced, 
and, moreover, the NLL approximation remains 
stable up to smaller scales than the fixed-order 
calculation. At the scale /i' « 1.3 GeV, which 
is close to the inverse Bohr radius, the NLL cor- 
rection vanishes. Furthermore, at /i" w 1.5 GeV, 
where a^^ = 0.319, the result becomes indepen- 
dent of fi; i.e., the NLL curve shows a local max- 
imum. This suggests a nice convergence of the 
logarithmic expansion despite the presence of the 
ultrasoft contribution with normalized at the 
rather low scale fj? /nib ~ 0.8 GeV. By taking 
the difference of the NLL and LL results at the 
local maxima as a conservative estimate of the er- 
ror due to uncalculated perturbative higher-order 
contributions, one gets E^hfs = 39 ± 8 MeV. A 
similar error estimate is obtained by the varia- 
tion of the normalization scale in the physically 
motivated soft region 1 — 3 GeV. 

For the case of charmonium, the use of pertur- 
bation theory is more doubtful, even though there 
has been some attempts in this direction recently 
[12] (see also [13] for the B^), yet one can try 
and see what comes out. This is specially impor- 
tant since even unqucnched attempts to obtain 
the HFS of charmonium from lattice undershot 
the experimental value by around 20% [14]. The 
results obtained in Ref. [4] are given in Fig. 2 
along with the experimental value 117.7±1.3 MeV 
[15]. The local maximum of the NLL curve corre- 
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Figure 1. HFS of IS bottomonium a.s a function 
of the renormalization scale /i in the LO (dotted 
line), NLO (dashed line), LL (dot-dashed line), 
and NLL (solid line) approximations. For the 
NLL result, the hand reflects the errors due to 
as{Mz) = 0.118 ±0.003. 



sponds to £;hfs = 104 MeV and a^^ = 0.534. We 
should emphasize the crucial role of the resumma- 
tion to bring the perturbative prediction closer to 
the experimental figure. Therefore, in Ref. [4] the 
whole difference of w 14 MeV between the pertur- 
bative prediction and the experimental value for 
the HFS of the ground state of charmonium was 
used to estimate the size of the non-perturbative 
effects. In any case, within the power counting as- 
sumed in Ref. [4], these non-perturbative effects 
are beyond the accuracy of this computation and 
are added to the errors. Taking into account that 
they are suppressed by the inverse heavy-quark 
mass as l/(asTOq)^ due to the multipole expan- 
sion, one obtains ~ 3.5 MeV for the typical size 
of the nonperturbative contribution to the HFS 
in bottomonium. For the estimate of the non- 
perturbative error, this number was multiplied 
by two, which was used above for the determina- 
tion of the theoretical error. These formulae has 
also been applied to n = 2 excited states. For 
bottomonium, one obtains i?hfs(2S')/i?hfs(l'S') = 
0.25. For charmonium, our perturbative estimate 
i?hfs(25')/£^hfs(lS') = 0.37 also reasonably agrees 
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Figure 2. HFS of IS charmonium as a function 
of the renormalization scale fi in the LO (dotted 
line), NLO (dashed line), LL (dot-dashed line), 
and NLL (solid line) approximations. For the 
NLL result, the band reflects the errors due to 
as{Mz) = 0.118 ± 0.003. The horizontal band 
gives the experimental value 117.7±1.3 MeV [15]. 



with the result 0.41 ± 0.03 of the recent experi- 
mental measurements [16]. Although one cannot 
rely on the (even NRG-improved) perturbative 
analysis of the excited charmonium states, the 
above agreement suggests that the nonperturba- 
tive effects are small, at least for the ground state. 

It has also been possible to give a good predic- 
tion for the HFS of the Bc{lS) system^ [5] 

M{B*)~M{Bc) =65 ±24 (th) t (fo,) MeV(2) 

Potential models appear to obtain slightly larger 
numbers [17]. As in the previous cases the in- 
clusion of logarithms appears to be a large ef- 
fect. In Fig. 3 from ref. [5], the HFS for the 
charm-bottom quarkonium ground state is plot- 
ted as a function of ly in the LO, NLO, LL, and 

■^In this case it is the pseudoscalar, Bc{lS), which has 
been measured, nevertheless not with enough precision to 
allow a good absolute determination of the B*(IS) mass. 
However, with improved statistics and precision of the Be 
data, this result can be considered as a prediction for the 
B* meson mass. 



NLL approximations for the hard matching scale 
value Vfi = 1.95 GeV. As we see, the LL curve 
shows a weaker scale dependence compared to 
the LO one. The scale dependence of the NLO 
and NLL expressions is further reduced, and, 
moreover, the NLL approximation remains sta- 
ble at the physically motivated scale of the in- 
verse Bohr radius, Cpotsmr ^ 0.9 GeV, where 
the fixed-order expansion breaks down. At the 
scale v' K, 0.85 GeV, which is close to the inverse 
Bohr radius, the NLL correction vanishes. Fur- 
thermore, at v" = 0.92 GeV, the result becomes 
independent of v; i.e., the NLL curve shows a lo- 
cal maximum corresponding to _Ehfs = 65 MeV, 
which is taken as the central value of the estimate. 
The NLL curve also shows an impressive stability 
with respect to the hard matching scale variation 
in the physical range rUc < vt < f^b and has a 
local maximum at Vh = 1.95 GeV, which is taken 
for the numerical estimates. All this suggests a 
nice convergence of the logarithmic expansion de- 
spite the presence of the ultrasoft contribution 
where as is normalized at the rather low scale 
l^h ^ 0.5 GeV. Let us discuss the accuracy of 
our result. For a first estimate of the error due to 
uncalculated higher-order contributions, we take 
9 MeV, the difference of the NLL and LL results 
at the local maxima. A different estimate can 
be obtained by varying the normalization scale in 
the physical range 0.8 < < 1.4 GeV. In this 
case the difference with the maximum is 16 MeV. 
Being conservative, we take this second number 
for our estimate of the perturbative error. Within 
the power counting assumed in Ref. [5], the non- 
perturbative effects are beyond the accuracy of 
the computation and added to the errors. They 
are also inferred using charmonium data in the 
same way than above and « 18 MeV for the er- 
ror due to the nonperturbative contribution to 
the HFS in Be was obtained. 

3. Production and Annihilation rates 

There has also been much progress on the re- 
summation of large logarithms appearing in the 
annihilation of a heavy quarkonium into leptons, 
photons or light hadrons, as well as its production 
in e+e^ or 77 collisions. 
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Figure 3. HFS for charm-bottom quarkonium as 
the function of the renormalization scale v in 
LO (dotted line), NLO (dashed line), LL (dot- 
dashed line ), and NLL ( solid line ) approximation 
for Vh = 1-95 GeV. For the NLL result the band 
reflects the errors due to as{Mz) = 0.118±0.003. 



The resummation of the large logarithms of the 
heavy quark velocity to all orders in as has been 
advocated as a tool to improve the behaviour 
of the perturbative expansion for tt threshold 
production [18]. Currently, the complete next- 
to-leading logarithmic (NLL) approximation for 
the production and annihilation rates is available 
[19,20]. In Ref. [21] a phenomenological analysis 
of the NLL result for the electromagnetic inclu- 
sive decays of the heavy quarkonium was made. 
The first attempt to go beyond the NLL approxi- 
mation [18] suggested a very good convergence of 
the logarithmic expansion. In particular, an ac- 
curacy of 2-3% was claimed for the cross section 
of tt threshold production. However, subsequent 
calculations of some next-to-next-to-leading log- 
arithmic (NNLL) terms [22], which had not been 
taken into account in Ref. [18], casted serious 
doubts on this estimate. Thus, the full calcula- 
tion of the NNLL corrections, which still remains 
elusive, is unavoidable to draw definite conclu- 
sions. In Ref. [23], it has recently been derived 



Figure 4. The spin ratio as the function of 
the renormalization scale v in LO=LL (dotted 
line), NLO (short-dashed line), NNLO (long- 
dashed line), NLL (dot-dashed line), and NNLL 
(solid line) approximation for the (would be) to- 
ponium ground state with Vh — rnt- For the 
NNLL result the band reflects the errors due to 
as{Mz) = 0.118 ± 0.003. Note that for the verti- 
cal axis the zero is suppressed. 



the complete NNLL result for the spin depen- 
dent part of the heavy quarkonium production 
annihilation rates, which includes the terms of 
the form In" for all n and applied to the 
heavy quarkonium phenomenology. In Figs. 4, 
5 and 6, the spin ratio is plotted as a function 
of 1/ in the various logarithmic and fixed-order 
approximations for the (hypothetical) toponium, 
bottomonium and charmonium ground states, re- 
spectively. As we see, in the second order the con- 
vergence and stability of the result with respect 
to the scale variation is substantially improved if 
one switches from the fixed-order to the logarith- 
mic expansion. We want to remark that the v 
dependence of the NLL approximation is slightly 
worse than at NLO. This is due to the artificially 
small I' dependence at NLO which is likely due 
to the fact that at this order only the hard scale 
enters. 
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Figure 5. The spin ratio as the function of the 
renormalization scale v in LO=LL (dotted line), 
NLO (short-dashed line), NNLO (long-dashed 
line), NLL (dot-dashed line), and NNLL (solid 
line) approximation for the hottomonium ground 
state with Vh = rub- For the NNLL result the hand 
reflects the errors due to as{Mz) = 0.118±0.003. 



Let us first consider the top quark case. In par- 
ticular, the ratio of the cross sections of the res- 
onance e^e~ — s- ti and 77 tt production. As 
one can see in Fig. 4, the logarithmic expansion 
shows perfect convergence and the NNLL correc- 
tion vanishes at the scale 1/ « 13 GeV, which is 
close to the physically motivated scale of the in- 
verse Bohr radius asmt/2. For illustration, at 
the scale of minimal sensitivity, ly = 20.2, GeV 
we have 



(Tres(e+e 



tt) 



,(77 ^ tt) 



^(1-0.132-0.018) . (3) 



However, it is not clear if the nice behaviour of 
the logarithmic expansion also holds for the spin- 
independent part of the threshold cross section. 
A possible problem is connected to the ultra- 
soft contribution, which is enhanced by the larger 
value of as at the ultrasoft scale. Whereas it is 
suppressed in the spin ratio by the fifth power 
of as, for the spin-independent part it already 



Figure 6. The spin ratio as the function of 
the renormalization scale v in LO=LL (dotted 
line), NLO (short-dashed line), NNLO (long- 
dashed line), NLL (dot-dashed line), and NNLL 
(solid line) approximation for the charmonium 
ground state with Vh = TOc- For the NNLL result 
the lower (yellow) hand reflects the errors due to 
as{Mz) = 0.118 ±0.003. The upper (green) hand 
represents the experimental error of the ratio [15] 
where the central value is given hy the horizontal 
solid line. 



contributes at 0{al) and can destabilize the ex- 
pansion. 

For bottomonium, the logarithmic expansion 
shows nice convergence and stability (c.f. Fig. 5) 
despite the presence of ultrasoft contributions 
with Os normalized at a rather low scale jra^. 
At the same time, the perturbative corrections 
are important and reduce the leading order re- 
sult by approximately 41%. For illustration, at 
the scale of minimal sensitivity, v — 1.295 GeV, 
we have the following series: 



r(T(15)-.e+e-) _ 1 

r(^,(15) 77) ~iQV 0-111) -(4) 

In contrast, the fixed-order expansion blows up 
at the scale of the inverse Bohr radius. 

So far we have discussed the perturbative cor- 
rections to the Coulomb-like quarkonium. How- 



ever, in contrast to the tt system, for bottomo- 
nium nonperturbative contributions can be im- 
portant. In our case the interaction of the quark- 
antiquark pair to the nonperturbative gluonic 
field is suppressed by v through the multipole ex- 
pansion in the same way than for the HFS com- 
putation. In this case however they only con- 
tribute in the N'^LL approximation, far beyond 
the precision of this computation. Note that the 
nonperturbative contribution to the decay rates 
ratio is suppressed by a factor in comparison 
to the binding energy and decay rates, where the 
leading nonperturbative effect is due to chromo- 
electric dipole interaction. Thus, by using the 
available experimental data on the T meson as 
input, one can predict the production and anni- 
hilation rates of the yet undiscovered rjh meson. 
In particular, one can predict the 77f,(lS') decay 
rate using the experimental value for the T(IS') 
decay rate and the following figure was obtained 
inRef. [23] 

r(?7f,(15) ^ 77) = 0.659 ±0.089(th.)ig;°i^(as) 
±0.015(exp.) keV, (5) 

where i' = 1.295 GcV, the scale of minimal sensi- 
tivity, was taken for the central value, the differ- 
ence between the NLL and NNLL result for the 
theoretical error and as{Mz) = 0.118 ± 0.003. 
The last error in Eq. (5) reflects the experimental 
error of r(T(15) -> e+e") = 1.314 ± 0.029 keV 
[15]. This value considerably exceeds the result 
for the absolute value of the decay width obtained 
in Ref. [21] on the basis of the full NLL analysis 
including the spin- independent part (see Fig. 7): 

r(?76(lS') ^ 77) = 0.35±0.1(th.)±0.05(as)KeV.(6) 

This can be a signal of slow convergence of the 
logarithmic expansion for the spin-independent 
contribution which is more sensitive to the dy- 
namics of the bound state and in particular to 
the ultrasoft contribution as it has been discussed 
above. On the other hand, the renormalon ef- 
fects [24] could produce some systematic errors in 
the pure perturbative evaluations of the produc- 
tion/annihilation rates. The problem is expected 
to be more severe for the charmonium case dis- 
cussed below. 
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Figure 7. Plot of r{7jbilS) 77) with LO/LL 
(dashed line), NLO (dot-dashed line) and NLL 
(dotted line) accuracy versus the renormalization 
scale V. 



We would also like to remark that the one-loop 
result for v = m\y overshoots the NNLL one by ap- 
proximately 30%. This casts some doubts on the 
accuracy of the existing as determination from 
the r(T light hadrons)/r(T e+e~) decay 
rates ratio, which gives as(jn\,) = 0.177 ± 0.01, 
well below the "world average" value [15]. The 
theoretical uncertainty in the analysis is esti- 
mated through the scale dependence of the one- 
loop result. Our analysis of the photon medi- 
ated annihilation rates indicates that the actual 
magnitude of the higher order corrections is most 
likely quite beyond such an estimate and the the- 
oretical uncertainty given in [15] should be in- 
creased by a factor of two. This brings the result 
for as into Icr distance from the "world average" 
value. 

For the charmonium, the NNLO approxima- 
tion becomes negative at an intermediate scale 
between agm-c and (c.f. Fig. 6) and the use of 
the NRG is mandatory to get a sensible pertur- 
bative approximation. The NNLL approximation 
has good stability against the scale variation but 
the logarithmic expansion does not converge well. 
This is the main factor that limits the theoretical 
accuracy since the nonperturbative contribution 
is expected to be under control. For illustration, 
at the scale of minimal sensitivity, v — 0.645 GeV, 



one obtains [23] 

r(J/^(15') ^ e+e- 
r(77,(15)^77) 



0.00001 

6 



0.513- 0.326) .(7) 8-10 



The central value is 2a below the experimental 
one. The discrepancy may be explained by the 
large higher order contributions. This should 
not be surprising because of the rather large 
value of as at the inverse Bohr radius of char- 
monium. For the charmonium HFS, however, 
the logarithmic expansion converges well and the 
prediction of the NRG is in perfect agreement 
with the experimental data. Thus one can try 
to improve the convergence of the series for the 
production/annihilation rates by accurately tak- 
ing into account the renormalon-related contribu- 
tions. One point to note is that with a potential 
model evaluation of the wave function correction, 
the sign of the NNLO term is reversed in the char- 
monium case [25] . At the same time the subtrac- 
tion of the pole mass renormalon from the per- 
turbative static potential makes explicit that the 
potential is steeper and closer to lattice and phe- 
nomenological potential models [26]. Therefore, 
the incorporation of higher order effects from the 
static potential may improve the agreement with 
experiment. Finally, we can not avoid mention 
that the NLL evaluation for the decay is able to 
reproduce the experimental value (see Fig. 8). 
Therefore, some extra work needs to be done to 
clarify these issues. 

4. Conclusions and outlook 

We have seen that the resummation of loga- 
rithms appear to have a large phenomenological 
impact in the heavy quarkonium physics. This is 
so in top-, bottom- or charmonium physics. 

First, the heavy-quarkonium HFS have been 
studied in the NLL approximation. The use of 
the NRG extends the range of /i where the per- 
turbative result is stable to the physical scale of 
the inverse Bohr radius. The resummation of log- 
arithms is found to be crucial to bring the pertur- 
bative prediction closer to the experimental figure 
of the HFS in charmonium (despite a priori un- 
suppressed nonpcrturbative effects), and to give 
reliable predictions for the rjbi^S) mass and the 
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Figure 8. Plot of T{rir{lS) 77) with LO/LL 
(dashed line), NLO (dot- dashed line) and NLL 
(solid line) accuracy versus the renormalization 
scale V. The horizontal line and its hand give the 
experimental value and its errors: r(?7c(l'5') 
77) = 7.2 ± 1.2 KeV [15]. 



Bc{lS) HFS. These results seem to indicate that 
the properties of the charmonium. Be and bot- 
tomonium ground states are dictated by pertur- 
bative dynamics. 

In the case of t-f production near threshold, the 
partial NNLL analysis made in Ref. [22] does not 
seem to show a very nice convergence (even if the 
absolute value of the corrections is small). Never- 
theless, being incomplete, such analysis is scheme 
dependent. In Ref. [23] a complete result with 
NNLL accuracy has been obtained for the ratio 
of the spin one and spin zero production. This is 
a physical result by itself and therefore scheme in- 
dependent. In this case a very nice convergence is 
found. Nevertheless, in this case the contribution 
due to the ultrasoft scale is suppressed. There- 
fore, it is premature to draw any definite conclu- 
sion for the convergence of the series. We are then 
eagerly waiting for the complete NNLL evalua- 
tion, which, even if difficult, is within reach. This 
is of utmost importance for the future determina- 
tions of the top mass and the Higgs-top coupling 
at the Next Linear Collider [27]. 

For the inclusive electromagnetic decays of the 
bottomonium and charmonium ground states the 
effects due to the resummation of logarithms ap- 



pear to be large and always improve the result 
compared with the finite order evaluations, yet 
the errors arc large and further work seems to 
be needed. Here as well, the complete NNLL 
evaluation would be of utmost importance to fur- 
ther clarify the physical picture. However, one 
should not forget the possible drawbacks of these 
determinations. There is an implicit dependence 
on the ultrasoft scale, which for the case of bot- 
tomonium and charmonium is quite low and is re- 
lated with non-perturbative effects, which should 
be eventually studied with the help of lattice or 
models. Moreover, renormalon effects could also 
play a role. 

Further work is ahead. No analysis has been 
made yet for the non-relativistic sum rules (to 
obtain them with NNLL accuracy would be a ma- 
jor step to obtain accurate determinations of the 
bottom and charm masses since they are strongly 
scale dependent), nor the impact of the resum- 
mation of logarithms in the determination of the 
bottom and charm masses from the ground state 
masses (which are known with NNLL accuracy 
[3,20]) estimated. 
Acknowledgments 

The author would like to acknowledge very pleas- 
ant collaborations with B.A. Kniclil, A. A. Pcnin, 
V.A. Smirnov and M. Steinhauser on which parts 
of the work reported here are based. 

REFERENCES 

1. A. Pineda and J. Soto, Nucl. Phys. B (Proc. 
Suppl.) 64, 428 (1998). 

2. A.H. Hoang, A.V. Manohar and LW. Stewart, 
Phys. Rev. D64, 014033 (2001). 

3. A. Pineda, Phys. Rev. D65, 074007 (2002). 

4. B.A. Kniehl, A.A. Penin, A. Pineda, V.A. 
Smirnov, and M. Steinhauser, Phys. Rev. 
Lett. 92, 242001 (2004). 

5. A.A. Penin, A. Pineda, V.A. Smirnov and M. 
Steinhauser, Phys. Lett. B593, 124 (2004). 

6. H. Stock, hep-ex/0310021. 

7. F. Maltoni and A.D. Polosa, hep-ph/0405082. 

8. M. A. Sanchis-Lozano, hcp-ph/0407320. 

9. N. Eicker et al, Phys. Rev. D57, 4080 (1998); 
T. Manke et al, Phys. Rev. D62, 114508 
(2000); X. Liao and T. Manke, Phys. Rev. 



D65, 074508 (2002). 

10. S. Godfrey and J.L. Rosner, Phys. Rev. 
D64, 074011 (2001); (E)65, 039901 (2002); 
F.J. Llanes-Estrada, S.R. Cotanch, A.P. Szc- 
zepaniak and E.S. Swanson, hep-ph/0402253. 

11. S. Narison, Phys. Lett. B387, 162 (1996). 

12. N. Brambilla, Y. Sumino, and A. Vairo, Phys. 
Lett. B513, 381 (2001); S. Recksiegel and Y. 
Sumino, Phys. Lett. B578, 369 (2004). 

13. N. Brambilla and A. Vairo, Phys. Rev. D62, 
094019 (2000). 

14. M. di Picrro et al., Nucl. Phys. Proc. Suppl. 
129, 340 (2004). 

15. S. Eidelman et al, Phys. Lett. B592, 1 
(2004). 

16. S.-K. Choi et al, Phys. Rev. Lett. 89, 102001 

(2002) ; 89, 129901(E) (2002); K. Abe et al, 
Phys. Rev. Lett. 89, 142001 (2002); J. Ernst 
et al, hep-ex/0306060; B. Aubert et al, hep- 
ex/0311038. 

17. E. Eichten and C. Quigg, Phys. Rev. D49, 
5845 (1994); E. Bagan, H. Dosch, P. Gos- 
dzinsky, S. Narison and J. Richard, Z. Phys. 
C64, 57 (1994); S. Gcrshtcin, V. Kisclcv, A. 
Likhoded and A. Tkabladze, Phys. Usp. 38, 
1 (1995); S. Godfrey, hep-ph/0406228. 

18. A.H. Hoang, A.V. Manohar, I.W. Stewart, 
and T. Teubner, Phys. Rev. Lett. 86 (2001) 
1951. 

19. A. Pineda, Phys. Rev. D66, 054022 (2002). 

20. A. H. Hoang and I. W. Stewart, Phys. Rev. 
D 67, 114020 (2003). 

21. A. Pineda, Acta Phys. Polon. B34 (2003) 
5295. 

22. A.H. Hoang, Phys. Rev. D 69 (2004) 034009. 

23. A. A. Penin, A. Pineda, V. A. Smirnov and 
M. Steinhauser, arXiv:hep-ph/0406175. 

24. E. Braaten and Y.Q. Chen, Phys. Rev. D57 
(1998) 4236; E-i&id. D59 (1999) 079901. 

25. A. Czarnecki and K. Melnikov, Phys. Lett. B 
519 (2001) 212. 

26. Y. Sumino, Phys. Rev. D65 (2002) 054003; 
S. Recksiegel and Y. Sumino, Phys. Rev. 
D65 (2002) 054018; A. Pineda, J. Phys. G29 

(2003) 371; T. Lee, Phys. Rev. D67 (2003) 
014020. 

27. M. Martinez and R. Miquel, Eur. Phys. J. C 
27 (2003) 49. 



